Short-time dynamics of a family of XY noncollinear magnets 
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Critical scaling and universality in the short-time dynamics for antiferromagnetic models on a 
three-dimensional stacked triangular lattice are investigated using Monte Carlo simulation. We have 
determined the critical point by searching for the best power law for the order parameter as a function 
of time and measured the critical exponents. Our results indicate that it is possible to distinguish 
weak first-order from second-order phase transitions and confirm that XY antiferromagnetic systems 
undergo a (weak) first order phase transition. 

PACS numbers: 75.10.Hk, 75.40.Cx, 75.40.Mg 



Frustrated systemaiS are characterized by competing 
interactions which may arise due to either disorder or 
geometry. The behavior of such systems is often unpre- 
dictable but the basic concepts of frustrated systems may 
provide insights into the physics of complex systems and 
have practical uses in areas ranging from microelectron- 
ics to drug dcliverji^iii^. Magnetic systems provide simple 
examples of frustration where exotic cooperative phases 
such as the 'spin glass', 'spin liquid' and 'spin ice' are 
found. During the past twenty five years a great deal of 
research effort has been put into investigating the nature 
of phase transitions in Heisenberg and XY frustrated sys- 
tems in three dimensions A^SiSiiSiii Particular attention 
has been devoted to Heisenberg and XY stacked trian- 
gular antiferromagnets that are commonly referred to as 
STA models. These models represent the simplest situa- 
tion of frustration induced by the geometry of the lattice 
leading to a critical behavior distinct from that encoun- 
tered in the usual ferromagnetic case. Indeed, in a tri- 
angular lattice, the competition due to antiferromagnetic 
interactions between nearest neighbour spins leads to a 
ground state with a planar spin configuration. In each 
elementary triangular cell the spins form a 120° structure 
with the vector sum of the three spins is equal to zero: 



Sr: =0 
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where the subscripts A, B, C label the sites at the corner 
of the elementary triangles shown in figure 1 . As a conse- 
quence the order parameter is no longer a simple vector 
but a matrix, a fact that has lead to the idea that these 
(noncollinear) frustrated magnets could belong to a new 
"chiral" universality classic. 

The XY STA model has been used to describe a 
great number of stacked triangular materials including 
CsCuCls, CsNiCla, CsMnIa and CsCuCls as well as the 
XY Helimagnets Ho, Dy and Tb. The experimental re- 
sults indicate that these materials exhibit second order 
phase transitions with the exception of CsCuCls where 
the transitioni^ is found to be weakly first order. The 
measured critical exponents exhibit scaling laws but vary 



from material to material which contradicts the basic 
idea of a unique set of critical exponents for all mate- 
rials described by the same model. In some experiments 
and also in some numerical simulations, the critical ex- 
ponent 77, also called the anomalous dimension, is nega- 
tive. This is forbidden if the theory which describes the 
transition is a unitary Landau Ginzburg Wilson (LGW) 
modeU4. Theoretical investigations using a perturbative 
renormalization group (RG) calculation up to high or- 
der predict the existence of a fixed point and, thus, the 
possibility of a second order phase transitio n^^i^^ . The 
varying critical exponents in this study are associated 
with a spiral-like RG fiow to a chiral, focus fixed pointi&. 
Non-perturbative RG methods (NPRG) predict a weak 
first order phase transition and attribute the appearance 
of scaling by a slowing down of the RG flow in the whole 
region of the coupling constant spac Q^i^° . 

The first numerical investigation of these STA mod- 
els using Monte Carlo methods indicated a second order 
phase transition with set of critical exponents possibly 
associated with a new chiral universality classic. Some 
subsequent numerical investigations have been performed 
on a modified version of the STA modelii, the STAR 




FIG. 1: Ground state of the STA system. The shaded trian- 
gles represent elementary plaquettes. 
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model, with the R representing a rigid constraint. In 
this model, the 120° structure of the ground state is lo- 
cally imposed at all temperatures. As a consequence the 
fluctuations of the spins within a triangular cell are sup- 
pressed while the fluctuations in the relative orientation 
of the disconnected triangular cells, or plaquettes, can 
still occur. Note that the STA and STAR models have 
the same symmetries and the "microscopic" changes per- 
formed are supposed to be irrelevant to the critical be- 
havior if it is universal. In fact, it was found that the 
STAR model exhibits a strong first order phase transi- 
tion thus raising doubts about the second order character 
of the phase transition occuring in all XY STA models. 
Finally, recent numerical studies of the STA model and 
its LGW formulation by ItakuraiAS, also indicate a first 
order phase transition for the STA XY model itself. 

In order to examine this effect of local rigidity in more 
detail, we introduced a generalized model in which we 
can continuously tune the local rigidity from the STA to 
the STAR hmitsi^i^a, 

^fW = -I]^^.S,.S,+r (Sa + Sb + Sc)' (2) 

(ij) plaquettes 

The interactions Jy are antiferromagnetic within the tri- 
angular layers and ferromagnetic between layers and have 
the same magnitude J = 1. The subscripts A, B, C label 
the three sublattices on the corners of each elementary 
triangle and the plaquettes refer to disconnected trian- 
gles as shown in figure 1. The parameter r imposes a 
constraint on the short wavelength fluctuations of the 
order parameter. Continuous changes in r from zero to 
infinity correspond to a continuous change from the STA 
to the STAR model. 

In our previous work with systems of linear sizes L < 
60 we found two different types of behavior: for r < 1.0 
the system exhibits a 'pseudo-critical' behavior whereas, 
for r > 1.0, a first order phase transition occurs. The 
critical exponents obtained in the r < 1.0 range appear 
to vary with the rigidity parameter r. This nonuniversal 
behavior is inconsistent with true critical behavior at a 
continuous phase transition for systems having the same 
symmetry of the order parameter. We concluded that the 
critical exponents are really 'pseudocritical' exponents 
and the observed scaling is 'pseudoscaling'. The esti- 
mated values of critical exponents are within the range of 
the experimentally observed critical exponents for ABX3 
compounds and Tb. In the range r > 1.0 we were able 
to estimate the value of latent heat for several values of 
r. We extrapolated the values of the latent heat to r = 
and we found a small but nonzero latent heat for the XY 
STA model which indicated a very weak first order phase 
transition. This behavior was confirmed by studying the 
energy probability distribution using much larger systems 
sizes L = 96, 138. Even larger sizes would be needed in 
the case of the Heisenberg modeU^. At negative values of 
r, the plaquettes are aligned ferromagnetically but inter- 
act antiferromagnetically and the symmetry of the order 
parameter is the same as at positive r. A special case 



occurs at r = —1/2 where we simply have a system of 
stacked Kagome layers. Additional degeneracies are ex- 
pected in this case. 

The standard equilibrium Monte Carlo approach re- 
quires extremely large lattice sizes L and long runs to 
properly sample statistically independent configurations. 
Hence, reaching a definite conclusion about the nature of 
the phase transition that occurs in STA models requires 
a different approach. This need is even more important 
for the Heisenberg STA which, from a numerical or theo- 
retical point of view, is expected to be closer to a second 
order behavior than the XY STA. For this reason we use 
an approach based on short time critical dynamicsSLS^. 
Universality and scaling have been observed in systems 
far from equilibrium. The initial state can be a high T 
state which is rapidly quenched to the critical tempera- 
ture or it can be a completely ordered state heated up 
to this temperature. Janssen et. aliSi showed that if a 
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FIG. 2: The z component of the order parameter mo as a func- 
tion of time for r = 4 and L = 72 with (a) an ordered initial 
state and {h) a disordered initial state for various tempera- 
tures. The solid line is obtained by quadratic interpolation 
and a least square fit to the expected power law behavior. 
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system is prepared at high T with an initial value of the 
order parameter tjiq and then quenched to the critical 
temperature, the time dependent order parameter ra{t) 
obeys the following scaling form after a small microscopic 
time 

m(t, T, L, mo) = &-'^/''m(^/&^ fo^/'^T, L/b, 6^«mo) (3) 

At the critical point (t = 0), this form predicts that 
m{t) increases as toq where 6 = (a;o — (3/v)/ z is a, new 
critical exponent and the other exponents /?, v, z are the 
usual critical exponents. Higher moments of the order 
parameter behave similarly. The critical temperature 
is located by starting from an initial value of toq a-t high 
T and then quenching to a lower T ^ Tc. A log-log plot 
of m(t) versus time should be linear and the slope yields 
6. An ordered state can also be used as an initial state to 
estimate Tc and in this case the slope yields f3/i^z . For 
a 2nd order transition these values of Tc should agree 
and the exponents should be universal. The derivative 
9t- lnm(t, T)|r=o should also exhibit power law behavior 
with the exponent l/vz. The advantage of this approach 
is that critical slowing down is not an issue since the 
measurements are carried out at short times. Averages 
are performed over different realizations of the initial val- 
ues of mo and thus time averages are replaced by sample 
averages and the full power of parallelized codes can be 
used. 

The ground state of the model has the spins on the 
corners of each elementary triangle arranged at 120° to 
one another. There are three sublattices as indicated in 
figure 1 and the order parameter can be chosen as 



mo = Sa - 0.5Si3 - 0.5Sc 



(4) 



We consider both a low temperature and high temper- 
ature initial state and then we follow the order param- 
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FIG. 4: The critical temperatures for both disordered and or- 
dered initial states plotted as a function of L for (a) r = 2 and 
(b) r — . The lines indicate the average Tc in each case and 
predict a ATc = 0.013(l)(r = 2) and ATc = 0.0043(l)(r = 0) 
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FIG. 3: The r — 4 critical temperatures for both disordered 
and ordered initial states plotted as a function of L. The 
lines indicate the average Tc in each case and predict a ATc — 
0.023(1) 



eter as a function of Monte Carlo time steps at some 
intermediate fixed temperature using the Metropolis al- 
gorithm. The results for each initial state are averaged 
over 10'^ — 20 X lO'^ trials depending on the linear size L. 

Figure 2 shows the behavior for a lattice of linear size 
L = 72 and r = 4 for several different temperatures start- 
ing from (a) the ordered state and (6) the disordered 
state. For the ordered state shown in figure 1 only one 
component of mo is non zero and we label this component 
as mo^. For the ordered initial state it has its maximum 
value and for all disordered initial states it has the ini- 
tial value 0.01. We follow the time dependence of moz as 
a function of time t for three temperatures and we use 
an interpolation scheme to determine the temperature 
Tc which yields the best power law behavior in the time 
range [100, 1000] as shown by the solid curve in the fig- 
ures. This intermediate time range was found to give the 
best power law behavior which only emerges after a time 
period which is long in microscopic termaSi. The error- 
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TABLE I: Results for the difFerence in critical temperatures 
ATc and the exponents for various values of the constraint 
parameter r. 



r 


AT, 


e 


P/vz 


l/vz 


P 





0.0043(1) 


0.081(4) 


0.218(1) 


0.79(2) 


0.27(1) 


2 


0.013(1) 


0.045(5) 


0.188(2) 


0.83(2) 


0.23(1) 


4 


0.023(1) 


0.028(5) 


0.169(1) 


0.82(1) 


0.20(1) 



bars on Tc are determined by the number of intermediate 
temperatures used in the interpolation scheme. We have 
performed the same calculations for sizes L = 60,72,90 
and the values of Tc are shown in figure 3 plotted ver- 
sus L. The values of Tc are independent of L in this 
time interval and the system displays hysteresis with 
the difference in the values of Tc corresponding to the 
ordered and disordered inital states yielding a value of 
ATc = 0.023(1). We have carried out the same proce- 
dure for smaller values of r. Figures 4(a), 4(b) show the 
behavior of Tc for r = 2 and r — respectively. In both 
cases the results indicate a weak first order transition and 
the values of the critical temperatures obtained here us- 
ing the present approach straddle those obtained using 
equilibrium methods^ii. Table I summarizes the results 
for the differences in Tc as determined from the ordered 
and disordered initial states as well as estimates for the 
various critical exponents obtained from the best power 
law dependence on t. The values of the exponents vary 
with the constraint parameter r which indicates nonuni- 
versal behavior. Using our measured values of (3/i>z and 
1 /vz we estimate the values of (3 given in the last column. 
The values of /? increase as r decreases in agreement with 
our previous study using equilibrium methods The 
value of /3 = 0.27(1) for r = is slightly larger than 
that predicted by previous numerical studies but is con- 



sistent with the value obtained in experiments on STA 
XY antiferromagnetsiS 

In summary, we have investigated the critical behavior 
of a family of XY noncoUincar magnets on the stacked 
triangular lattice geometry using the short-time dynam- 
ics approach. The critical temperatures obtained using 
this approach straddle those obtained previously using 
equilibrium methods and indicate that the transition is 
accompanied by hysteresis. The critical exponents are 
found to vary with the constraint parameter r. Since this 
parameter does not change the symmetry of the model 
the exponents are non-universal. Our results strongly 
suggest that the phase transition of STA XY antifer- 
romagnets is weakly first order in agreement with the 
NPRG field theory predictions and with our previous 
equilibrium Monte Carlo results. The method used here 
has the advantage that scaling behavior emerges at rel- 
atively short times and also for smaller sizes since our 
values of Tc arc almost independent of L. The results in- 
dicate that, for the STA XY materials, experiments need 
to carried out at reduced temperatures r <C 10~^ in or- 
der to identify the true weak first order nature of the 
transition. The present method also provides a way to 
study the question of a even weaker first order transition 
for STA Heisenberg materials. 
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